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Abstract 


Hunted to near extinction in the late 19th century, the endangered and endemic Hawaiian monk seal (Monachus schauinslandi) 
exhibits low variation at all molecular markers tested to date. Here we confirm extreme paucity of genetic diversity, finding 
polymorphisms at only 8 of 154 microsatellite loci tested (143 novel species-specific loci, 10 loci from Antarctic seals, and 
1 previously characterized locus). This screening revealed unprecedentedly low levels of allelic diversity and heterozygosity 
(A = 1.1, He = 0.026). Subsequent analyses of 2409 Hawaiian monk seals at the 8 polymorphic loci provide evidence for 
a bottleneck (P = 0.002), but simulations indicate low genetic diversity (H. < 0.09) prior to recorded human influence. 
There is little indication of contemporary inbreeding (Fg = 0.018) or population structure (K = 1 population). Minimal 
genetic variation did not prevent partial recovery by the late 1950s and may not be driving the current population decline 
to ~1200 seals. Nonetheless, genotyping nearly every individual living during the past 25 years sets a new benchmark for 


low genetic diversity in an endangered species. 


Genetic variation is the raw material of evolution and 
a prerequisite for adaptation via natural selection (Fisher 
1930). Measures of genetic diversity are often used to 
estimate individual fitness and the potential for population 
persistence (Coltman et al. 1999; Hansson and Westerberg 
2002; Reed and Frankham 2003). Genetically depauperate 
species may have a reduced ability to mount an effective 
defense against pathogens (Hawley et al. 2005) or to adapt 
to environmental changes (Hoffmann and Parsons 1997), 
thereby increasing the risk of extinction (Mills and Smouse 
1994; Lacy 1997; Frankham et al. 1999, 2002). Although 
extinction may be a demographic process (Lande 1988; 
Allendorf and Luikart 2007), meta-analysis of laboratory and 
wild case studies indicates that inbreeding depression and 
low genetic diversity cannot be ignored (Frankham 2005) as 
most species are not driven to extinction prior to being 
impacted by genetic factors (Spielman et al. 2004). 
Measures of genetic diversity, such as proportion of 
polymorphic microsatellite loci (P), number of alleles per 
locus or allelic diversity (4), and mean percent heterozy- 
gosity (H), reflect previous demographic events. In a meta- 
analysis of 108 mammalian species, Garner et al. (2005) do 
not find a consistent trend in the genetic diversity of healthy 


populations, but there is a pervasive and consistent 
reduction of genetic variability in populations that have 
experienced a rapid and severe decline in population size. 
However, such bottlenecks do not always produce the 
predicted genetic signature (Busch et al. 2007) and are not 
the only cause of low genetic diversity. Amos and Harwood 
(1998) assert that inbreeding, population structure, hetero- 
zygote mutation bias, and genome-wide selective sweeps 
also contribute alone or in synergy to the low variability 
detected in many endangered species. 

The northern elephant seal (Mirounga angustirostris) was 
the first documented case of low genetic diversity in 
response to near extinction. In a seminal paper, Bonnell and 
Selander (1974) report unprecedented homogeneity at 
24 allozyme loci across 159 individuals. Hoelzel et al. (1993) 
confirm these results with another 43 allozyme loci and also 
find low variation at mitochondrial loci, with 2 haplotypes in 
40 individuals. At exon 2 of the class II major histocom- 
patibility complex DQB locus, expected to be highly 
polymorphic in mammals, Hoelzel et al. (1999) found only 
2 alleles in 69 individuals. Finally, in a survey of 39 
microsatellite loci, Garza (1998) identifies 14 polymorphic 
loci with 2—4 alleles each. A 19th century population 
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bottleneck, the result of overhunting, may account for much 
of the homogeneity: a simulation model based on extensive 
demographic data indicates a bottleneck of less than 30 seals 
and 20-year duration (Hoelzel et al. 1993). High reproductive 
skew, low reproductive rates, and demographic stochasticity 
also may have contributed to the low genetic variation 
(Hoelzel 1999). 

The Hawaiian monk seal (Monachus schauinslandi) endured 
a similar history of overexploitation. Originally occurring 
throughout the Hawaiian archipelago (Figure 1), it was likely 
extirpated from the main Hawaiian Islands by Polynesian 
colonizers 1500-1600 years ago (Bellwood 1978; Baker and 
Johanos 2004). On arrival of the first European sailors in 
the 19th century, the Hawaiian monk seal was hunted to 
neat extinction at the 6 primary Northwestern Hawaiian 
Islands subpopulations (French Frigate Shoals, Laysan, 
Lisianski, Pearl and Hermes Reef, Midway Atoll, and Kure; 
Ragen 1999). Although historical counts of total population 
size ate not available, records indicate an abundance of seals 
up to the year 1857 (Hiruki and Ragen 1992), no or few seals 
at most islands by 1893 (Ragen 1999), and a “large number” 
at Kur and Pearl and Hermes Reef by 1915 (Hiruki and 
Ragen 1992). In 1958, mean counts of seals on the beach 
(an index of abundance) had recovered to 916 individuals, 
1 year of age or older (nonpups; Rice 1960). The mean 
count has since declined to 293 nonpups, and the most 
recent estimate of total abundance is 1247 seals (Carretta 
et al. 2007). Extensive monitoring of the Hawaiian monk 
seal over the past 25 years has yielded over 2400 genetically 
sampled individuals, which represents a large proportion of 
the entire species. 

Mitochondrial sequencing and multilocus fingerprint 
analyses indicate low genetic diversity in the Hawaiian 
monk seal (Kretzmann et al. 1997), and Aldridge et al. 
(2006) report exceptional uniformity in major histocompat- 
ibility class I genes. Of 18 pinniped species tested, the 
Hawaiian monk seal exhibits the lowest level of genetic 
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diversity when assessed across 20 microsatellite loci isolated 
from gray (Halichoerus grypus, harbour (Phoca vitulina), 
southern elephant (Mirounga leonina), and South American 
fur (Arctocephalus australi) seals (Gemmell et al. 1997). 
Among these 20 loci, 3 are polymorphic in the Hawaiian 
monk seal (Gemmell et al. 1997), but only 1 locus (Hg6.3; 
Allen et al. 1995) conforms to expectations of Mendelian 
inheritance and can be used in population-level analyses 
(Kretzmann et al. 2001). 

Low microsatellite diversity is often attributed to 
ascertainment bias, that is, loci isolated in one species are 
comparatively invariant in distantly related species (Ellegren 
et al. 1995). Reviewing genetic diversity of 108 mammalian 
species, cross-species amplification results in consistently 
and significantly reduced mean heterozygosity and number 
of alleles per locus (Garner et al. 2005). The success of 
cross-species microsatellite amplification is generally related 
to the evolutionary distance from the focal species; for 
species diverging 10-20 million years ago, approximately 
25% of primer sets will amplify polymorphic loci (Gemmell 
et al. 1997; Primmer et al. 1996). Fossil evidence suggests 
that monachines (southern seals including the monk, 
elephant, and Antarctic seals) split from phocines (the 
northern seals including the harbour and gray seals) 15—17 
million years ago (Demere et al. 2003; de Muizon 1982). 
Molecular data indicate a divergence of the monk seals and 
elephant seals 11-16 million years ago (Fyler et al. 2005; 
Arnason et al. 2006). Therefore, one would expect approx- 
imately 5 of the loci tested above to be polymorphic in monk 
seals due to ascertainment bias alone, and Gemmell et al. 
(1997) found 3 to be variable in the Hawaiian monk seal. 

The paucity of variable loci has prevented use of genetic 
data in management initiatives. Population structure remains 
equivocal, with results based primarily on one microsatellite 
locus (Kretzmann et al. 2001). Identifying individuals 
genetically, confirming maternity, and elucidating paternity 
would provide measures of male and female reproductive 
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Figure I. Range of Monachus schauinslandi in the Hawaiian archipelago. 
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rates, which are priorities to wildlife managers. Finally, it is 
currently unknown whether inbreeding depression contrib- 
utes to decline of the species. To address these questions, 
additional variable loci are needed. 

Here we further investigate the genetic diversity of the 
Hawaiian monk seal. We first test for cross-species 
amplification using 10 microsatellite loci isolated from 
Antarctic seals (Davis et al. 2002). We then isolate and test 
143 new microsatellite loci from the Hawaiian monk seal and 
genotype 2409 individuals at 7 novel and 1 previously 
published (Hg6.3) polymorphic loci. We evaluate possible 
causes of low genetic diversity, including population bottle- 
neck, inbreeding, and population subdivision. Building on 
previous studies of nuclear and mitochondrial DNA diversity 
(Gemmell et al. 1997; Kretzmann et al. 1997, 2001), our survey 
of 154 loci provides a compelling portrait of depleted genetic 
diversity across the genome of the Hawaiian monk seal. 


Materials and Methods 


Tissue samples from Hawaiian monk seals (NV = 2409) were 
collected opportunistically and during annual population 
assessments in the main Hawaiian Islands (2000-2007; N = 
54) and at 7 sites in the Northwestern Hawaiian Islands 
(1980-2007): Nihoa (7), French Frigate Shoals (N = 766), 
Laysan (656), Lisianski (310), Pearl and Hermes Reef (260), 
Midway (134), and Kure (222). Seals were tagged by 
punching a hole in their flippers; the resulting tissue plugs 
were preserved in 20% dimethyl sulfoxide salt-saturated 
solution, 95% ethanol, or frozen in liquid nitrogen 
(Henderson and Johanos 1988). Genomic DNA was 
extracted from all samples using the DNeasy Tissue Kit 
(Qiagen, Valencia, CA). 

To isolate microsatellite loci from the Hawaiian monk 
seal, we followed the protocol of Glenn and Schable (2005) 
with modifications. Genomic DNA was extracted from 
captive seal blood using the DNeasy Tissue Kit. Approx- 
imately 5 ug of high—molecular weight DNA was digested 
using the blunt-end cutting enzyme 4%] to avoid methylated 
CG sequences common to mammals. Resulting fragments 
containing microsatellite motifs were hybridized to biotiny- 
lated oligonucleotides and captured using magnetic beads 
(Dynal, Oslo, Norway). A T-tailed vector was created by first 
digesting a pZErO-2 vector (Invitrogen, Carlsbad, CA) with 
the enzyme EcoRV (Promega, Madison, WT) and then adding 
a T-tail by performing a 2-h extension (72 °C) with 100 mM 
2'deoxythymidine 5’-triphosphate and polymerase (Bioline, 
San Clemente, CA). The fragments were amplified, ligated 
into the T-tailed vector, and used to transform o-Select Gold 
Efficiency Competent Cells (Bioline, Randolph, MA). 
Recombinant clones were selected at random, amplified, 
and sequenced on an ABI 3730XL DNA Analyzer (Applied 
Biosystems, Foster City, CA) by Macrogen, Inc. (Seoul, South 
Korea). Of 624 positive clones, 470 contained microsatellite 
sequences. The software PRIMER 3 (Rozen and Skaletsky 
2000) was used to develop primers for 143 putative loci that 
exceeded 10 short tandem repeats. To gain insight into the 
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genomic coverage and possible identity of these 143 novel 
loci, sequences were blasted to the dog genome in GenBank 
(http://www.ncbi.nlm.nih.gov/). 

To test for polymorphism, we analyzed 8 Hawaiian monk 
seal DNA samples (1 from each Northwestern Hawaiian 
Islands subpopulation and 2 from the main Hawaiian Islands) 
at 143 novel microsatellite loci and a subset of microsatellite 
loci developed from 3 species of Antarctic seals (Weddell, 
Leptonychotes weddellii, leopard, Hydrurga leptonyx; and crabeater, 
Lobodon carcinophagus): Lw8, Lw10, Lw11, Lw16, Lw18, H14 
H115, H116, H120, and Lc5 (Davis et al. 2002). The 
polymerase chain reaction (PCR) protocol consisted of 15 
min of initial denaturation at 95 °C, followed by 35 cycles of 
denaturation (94 °C, 30 s), annealing (56—60 °C, 30 s), and 
extension (72 °C, 30 s), with a final extension (72 °C, 30 min) 
of the following 10-ul reaction: 2x Biomix Red (Bioline), 
5 uM each primer, and 30-50 ng genomic DNA. We ran 5 pl 
of amplified product on a 3% agarose gel, stained with 
GelStar (Lonza, Basel, Switzerland) and visualized the bands 
using a Molecular Image Gel Documentation XR System 
(Bio-Rad, Hercules, CA). 

Primers amplifying the 7 novel polymorphic loci and 
1 previously isolated locus (Hg6.3; Allen et al. 1995) were 
fluorescently labeled and combined in a single multiplex 
PCR using the following 6-ul reaction: 1x PCR mix 
(Qiagen), 2 UM primer mix, and 30-50 ng genomic DNA. 
The PCR protocol consisted of 15 min initial denaturation 
at 95 °C, followed by 35 cycles of denaturation (94 °C, 30 s), 
annealing (60 °C, 30 s), and extension (72 °C, 30 s), with a 
final extension (60 °C, 30min) to ensure the addition of a 
terminal adenine. Amplified products were tun on an 
ABI 3100 Genetic Analyzer (Applied Biosystems) by the 
Experimental Program to Stimulate Competitive Research 
Sequencing Facility at the Hawaii Institute of Marine 
Biology and scored using the software GENEMAPPER 4.0 
(Applied Biosystems). 

We scored all products 2 or more times to ensure 
accuracy. Global tests for Hardy-Weinberg equilibrium and 
linkage disequilibrium were conducted using FSTAT 2.9.3.2 
(Goudet 1995). Mendelian inheritance of the alleles was 
assessed using known mother—pup relationships. PCR error 
as a result of null alleles and large allele dropout was tested 
in MICROCHECKER 2.2.3 (Van Oosterhout et al. 2004). 
Scoring and human error was estimated by duplicating 
analyses (from extraction to scoring) for 72 randomly 
chosen individuals. 

The average number of alleles per locus (4) and 
unbiased heterozygosity (H; Nei 1987) were calculated 
using EXCEL MICROSATELLITE TOOLKIT 3.1 (Park 
2001). We also used the program to estimate polymorphism 
information content (PIC; the probability that the parental 
genotype may be deduced from rare alleles found in its 
progeny) and evaluated each locus as highly informative 
(PIC > 0.5), reasonably informative (0.5 > PIC > 0.25), or 
slightly informative (PIC < 0.25; Botstein et al. 1980). 

We assessed the impact of the 19th century population 
bottleneck on genetic diversity using the software 
BOTTLENECK 1.2.02 (Piry et al. 1999). In a bottlenecked 
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population, rare alleles are the first to be lost, lowering the 
mean number of alleles per locus. Heterozygosity is less 
affected, producing a transient excess in heterozygosity 
relative to that expected given the resulting number of 
alleles (Cornuet and Luikart 1996). Because we scored less 
than 20 microsatellite loci, as required by the strict stepwise 
mutation model (SMM), we performed 10 000 iterations of 
the Wilcoxon signed-rank test with the infinite allele model 
(SMM; Kimura and Crow 1964) and two-phase mutation 
(TPM) model (TPM incorporated 70% stepwise and 30% 
multistep mutations; DiRienzo et al. 1994). We also applied 
the M ratio method, which measures the number of alleles 
relative to the overall range in allele size using the software 
MPVAL (Garza and Williamson 2001). We ran simulation 
tests (10 000 replicates in CRITICALM) to provide average 
and lower 95% confidence intervals for M ratio values, given 
0 = 4 and default values for a TPM microsatellite mutation 
model. 

Generally, a bottleneck must persist over several 
generations to impact heterozygosity significantly. Following 
Amos and Balmford (2001), we considered various 
scenarios to describe the loss of heterozygosity relative to 
time and effective population size using the following 
formula: 


H, = (1 = 1/2N,)’, 


where Họ is the initial heterozygosity and H, is the 
heterozygosity ¢ generations after a decline to size Ne 
(James 1970). H, was calculated from the observed 
genotypes. We estimated Ho over several generations (¢ = 
1—4 generations or 15—60 years) and for various effective 
population sizes (Ne = 2-50). Due to the lack of historic 
total abundance information, we calculated a minimum 
postbottleneck population size from historic beach counts 
and records of monk seals available from the early 20th 
century (Hiruki and Ragen 1992), assuming that following the 
bottleneck the population grew at the maximum growth rate 
observed for the species (7% per year; Carretta et al. 2007). 

Contemporary inbreeding (calculated as Frs, a measure of 
within population heterozygote deficit, with 95% confidence 
intervals by bootstrapping over all loci) was evaluated using 
FSTAT. Population subdivision was assessed in STRUC- 
TURE 2.0 (Pritchard et al. 2000), which assigns individuals 
into K populations to achieve Hardy-Weinberg and linkage 
equilibrium. We used the default settings (admixture model, 
correlated allele frequencies, and A = 1) with a 10 000 burn- 
in length and 100 000 simulations to test K = 1—10 with 
3 repetitions each, verifying convergence of alpha and 
likelihood values. The program POWSIM 4.0 (Ryman and 
Palm 2006) was used to evaluate the statistical power 
provided by the 8 polymorphic loci to detect population 
structure in the Hawaiian monk seal. 


Results 


Only 7 of the 143 microsatellite loci isolated from the 
Hawaiian monk seal were polymorphic in the initial survey 
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(GenBank; accession numbers EU913763—EU913769). We 
followed a standard protocol (Glenn and Schable 2005) with 
modifications that were unlikely to reduce variability. 
Though resolution is limited using 3% agarose gels to 
assess polymorphism, we were able to detect alleles differing 
by 2 (locus Ms647) and 4 (loci Ms15 and Ms663) base pairs, 
and subsequent scoring of 2401 individuals on an automated 
sequencer revealed only 2 additional alleles not observed in 
the initial screening. All Antarctic seal microsatellite primer 
pairs yielded amplification products in the Hawaiian monk 
seal, but none was polymorphic. Of 143 microsatellite loci 
isolated from the Hawaiian monk seal, 47 could not be 
matched to sequences within the dog genome, whereas 96 
were aligned to sequences on dog chromosomes X, 1—9, 11, 
13, 15-18, 20-27, 29-31, and 33-38. No further analysis 
was performed on the monomorphic markers. 

A total of 8 polymorphic loci (7 novel and previously 
published Hg6.3) were used to genotype 2409 individuals 
from the main and Northwestern Hawaiian Islands (Table 1). 
No loci showed significant deviation from Hardy—Weinberg 
expectations. Due to excessive genome-wide homozygosity, 
polymorphism at 8 loci may have been maintained by direct 
or indirect balancing selection. In all, 5 of the 8 polymorphic 
loci could not be matched to sequences in the dog genome, 
and the remaining 3, Hg6.3, Ms504, and Ms663, were 
matched to sequences with no known gene product. Two 
pairs of loci (Hg6.3 and Ms23, P = 0.027; Ms504 and Ms23, 
P = 0.033) exhibited linkage disequilibrium, which was not 
significant after Bonferroni correction. In the 23 mother— 
pup comparisons, all shared at least one allele for all loci, 
supporting Mendelian inheritance. There was no evidence 
for null alleles, large allele dropout, or scoring error. 
Seventy-two specimen pairs with identical identification 
numbers (tissue samples collected independently from the 
same seal) were independently extracted and genotyped, to 
test for human error in analysis. Two specimen pairs 
displayed discrepancies: each contained one mismatched 
allele at a single locus, an error rate of 0.2% per allele. 

Genetic diversity was extremely low. Evaluated across all 
loci, the average number of alleles per locus was A = 1.1 
(+0.79) and the overall unbiased heterozygosity was H. = 
0.026 (£0.12). Excluding 146 monomorphic loci, the 
number of alleles per locus varied from 2 to 9, with a mean 
of A = 3.5 (£2.62) and an overall unbiased heterozygosity 
H. = 0.49 (+£0.061). PIC values indicated 2 highly 
informative (PIC > 0.5) and 6 moderately informative 
(0.5 > PIC > 0.25) loci (Table 1). 

Analyzing the data following 2 models of microsatellite 
mutation (SMM and TPM), tests for heterozygosity excess 
were significant (P = 0.002), indicating a population 
bottleneck (Table 1). We acknowledge that the assumptions 
of these analyses are not met, namely our population is likely 
not at mutation-drift equilibrium and our loci do not 
conform to either model of microsatellite mutation. 
Therefore, we considered another bottleneck test. Simulation 
tests of the M ratio provided an average (M = 0.84) and lower 
95% confidence interval (M, = 0.74), describing a critical 
value above which a bottleneck is indicated. Our data indicate 
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0.008 
0.043 
0.021 
-0.020 
0.043 
0.029 
0.015 
0.010 
0.018 


Fis 


M ratio 

00 

.00 
0.60 

00 
0.20 

.00 
1.00 
0.33 
0.76 


0.002 


86/1.48 
38/0.94 
26/111 
54/141 


(He — Heq)/SD IAM/TPM 


1.36/0.88 
2.40/2.16 
1.10/0.98 
2.22/2.00 


P 


No. of alleles 


2 
3.5 


0.62 
0.74 
0.69 
0.26 
0.30 
0.46 


PIC 


0.67 
0.55 
0.75 
0.47 
0.29 
0.47 
0.29 
0.35 
0.48 


Ho 


He 
0.68 
56 
0.78 
0.49 
0.29 
0.46 
0.31 
0.36 
0.49 


0. 


Size (bp) 
340-370 
158 162 
203-315 
297-317 
308 326 
115 117 
290 294 
227 237 


TTAGTGTGGAGTCACTTAGG 


CGC 


F: TCAACTTCTCAA’ ‘AGGATTCACA 


F: CCAAAGCCTATTTCTTTCAATCC 
R: AGCAGAGGCCCTAAGACAGG 
F: CTGAATTCATGCTGTATCTTGG 
R: GCAAAAAGGGATGAGCCATA 
F: CAGGGGACCTGAGTGCTTATG 
R: GACCCAGCATCAGAACTCAAG 


R: GIGCTTGGGACATGATGG 


F: 
R: GTCATTCCCTAGTGGTAAAGACTC 


R: GTGAGATGAATGCCCTTTGG 
F: GACTGGTAATTTACGCCCTACC 
R: AAGTGTTGGGTTGAAAATTGG 
F: ATCAGCTATCAGGGGTAGGG 
F: GAACTCCAAACAGCCATTCC 

R: CCTGCTCCTTCTTTCTGATCC 


Microsatellite loci isolated from the Hawaiian monk seal, Monachus schauinslandi 
Primer sequence (5’-3’) 


(CCTT)6CCCT(CCTT)6 
(GAAA)9GGAA(GAAA)8 


Repeat motif 
(GAAA)18 
(Allen et al. 1995) 


(TG)14 
(TC) 


(GT)13 
(AAG) 24 


F, forward; R, reverse; IAM, infinite allele model. Locus name is an abbreviation of the genus and species; the number denotes the clone ftom which the primers were developed. Repeat motif and primer sequence information are 
reported for a single individual. PCR product size represents the range of allele sizes. Unbiased expected and observed heterozygosity (H. and Ho; Nei 1987), PIC (Botstein et al. 1980), and number of alleles (overall statistics in bold) are 
based on 2409 individuals. Excess heterozygosity values are the difference between expected and actual heterozygosity divided by the standard deviation, for 2 models of microsatellite mutation. M ratio values and inbreeding coefficients 


(Fis) are measured at each locus and overall (in bold). 


Table I. 
Locus 
Ms15 
Ms23 
Ms265 
Ms504 
Ms647 
Ms663 
Hg6.3 
Overall 


Ms9 
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a bottleneck (M = 0.76); however, the signal was weak (less 
than the simulation average), which may indicate rapid 
population recovery (Garza and Williamson 2001). 

We estimated the minimum (A) = 0.030) and maximum 
(Ho = 0.094) initial heterozygosity after a bottleneck of 
different durations (£ = 1—4 generations) and effective 
population sizes (NV. = 2-50). After a single-generation 
bottleneck (~15 years) in which population decline was 
minimal (Ne = 50), there is essentially no loss of 
heterozygosity (Ho = 0.030). An extreme bottleneck (NV. = 
2) for 4 generations was required to reduce heterozygosity 
by 67% (from Ho = 0.09 to H, = 0.03). 

The inbreeding coefficient was low over all polymorphic 
loci (Fis = 0.018; 95% CI —0.001 to 0.023) and at each 
locus (Table 1). We found evidence for a single population 
cluster in STRUCTURE. For K = 1, estimates of posterior 
probability approached one, whereas for K = 2-10, the 
values approached zero. For K = 2-10, there was high 
variance among A values, and all individuals were roughly 
symmetrical and admixed over K populations (i.e., no strong 
population assignments), consistent with a lack of popula- 
tion structure. It is unlikely that such results reflect low 
resolving power of our markers as a simulation test 
(POWSIM; Ryman and Palm 2006) revealed a 100% 
probability of detecting divergence among Hawaiian monk 
seal populations if the true amount of differentiation 
corresponds to an gy value of 0.01 or greater. 


Discussion 


The Hawaiian monk seal (M. schauinslandi) exhibits extremely 
low genetic diversity. Analyzing 154 microsatellite loci, we 
find unprecedented levels of allelic diversity (4 = 1.1) and 
heterozygosity (H. = 0.026), which in a species with 17 
chromosomal pairs, can be considered a relatively thorough 
sampling of the genome. Only 7 of the 143 novel 
microsatellites isolated in this study are polymorphic. Four 
of these, plus the gray seal locus Hg6.3 (Allen et al. 1995), are 
diallelic, and the average heterozygosity of the 8 variable loci 
is low (H. = 0.49 + 0.061; Table 1). 

Our survey confirms the initial reports of low genetic 
diversity in the Hawaiian monk seal at 20 polymorphic 
phocid loci (Gemmell et al. 1997) and an additional 7 loci 
isolated from the harbour seal, found to be monomorphic in 
the Hawaiian monk seal (Kretzmann et al. 2001). Our low 
allelic diversity (4 = 1.1) is similar to the previous studies 
(A = 1.0-1.3; Table 2). We find slightly higher poly- 
morphism than Kretzmann et al. (2001) but lower than 
Gemmell et al. (1997). Note that these studies only 
considered loci previously found to be polymorphic in at 
least one other seal species, whereas we include 143 novel 
microsatellite loci, developed specifically for the Hawaiian 
monk seal. When we consider only polymorphic loci, our 
average heterozygosity (H. = 0.49) falls within the range of 
previously reported observed (H, = 0.1—0.8) and expected 
(H. = 0.1-0.53) heterozygosities at 6 Northwestern 
Hawaiian Islands subpopulations (Kretzmann et al. 2001). 
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Table 2. Comparison of measures of microsatellite diversity, including proportion of loci which ate polymorphic (P), average 


heterozygosity (H. which is not significantly different from M, in 


examples listed), and average number of alleles per locus (4) 


Sample Estimated 

P He A size census size Reference 
Hawaiian monk seal Monachus schauinslandi 8/154 0.03 1.1 2409 1247 This study 

3/20 — 13 5 Gemmell et al. (1997) 

0/7 — 10 5 Kretzmann et al. (2001) 
Mediterranean monk seal Monachus monachus 15/39 0.16 1.5 41-52 300 Pastor et al. (2004) 
Northern elephant seal Mirounga angustirostris 14/41 0.13 1.5 80—160 127 000 Garza (1998) 
Cheetah Acinonyx jubatus 7/10 0.39 3.4 10 15 000 Menotti-Raymond and O’Brien (1995) 
North Atlantic right whale Eubalaena glacialis 117/27 0.32 3:2 209 300 Waldick et al. (2002) 
Northern hairy-nosed wombat Lasiorhinus krefftii 9/16 0.27 1.8 28 70 Taylor et al. (1994) 
Ethiopian wolf Canis simensis 6/9 0.21 2.4 20 500 Gottelli et al. (1994) 


Taxa represent the most genetically depauperate species and the lowest values reported for each species based on tables from Pastor et al. (2004) and 


Frankham et al. (2002). Some studies could not be included in the table because key values, such as number of monomorphic loci or global H., were not 


reported. 


Therefore, we conclude that previously reported low genetic 
diversity is not a result of ascertainment bias but rather 
characteristic of the Hawaiian monk seal. 

The extremely low diversity reflects a population 
bottleneck. Using heterozygosity excess and M ratio 
analyses, we find molecular evidence for population 
depletion, corroborating historical accounts of intense 
hunting of the Hawaiian monk seal during the 19th century 
(Ragen 1999). 

Records from European exploration of the Northwest- 
ern Hawaiian Islands document seals as early as 1805. There 
are reports of at least 4 shipwrecks from 1842 to 1870 in 
which the crew consumed seals (V = 10—60-+) for survival. 
A vessel visiting all islands in 1859 returned to Honolulu 
with 1500 skins and 240 barrels of seal oil (Brooks 1860). By 
the late 1880s and early 1890s, at least 2 colonies appear to 
have been extirpated (Laysan and Midway) and the others 
severely depleted (Ragen 1999). Though incomplete, records 
suggest that the early 1890s marked the nadir of the 
bottleneck. 

Using these accounts as a guide, we infer prebottleneck 
heterozygosity (Ho = 0.03—0.09) by assessing severe to weak 
population decline (NV. = 2-50) persisting for a range of 
durations (1—4 generations). At the most extreme bottleneck 
(N: = 2) for a duration of 4 generations (~60 years), there 
would have been a ~67% reduction in heterozygosity, 
which still corresponds to surprisingly low levels prior to 


exploitation (Họ = 0.09). Furthermore, this estimate is 
almost certainly an overestimate of heterozygosity because it 
is unlikely that M, = 2 for 60 years. Based on accounts 
compiled by Hiruki and Ragen (1992), the species had 
recovered to at least 184 individuals by the early 1920s 
(Table 3). Projecting back in time and assuming that the seal 
population grew at 7% per year (the maximum rate 
observed for the species), there would have been at least 
23 seals that survived in 1893, the last year when 
a considerable number of seals (60—70 at Laysan Island) 
were teportedly killed (Hiruki and Ragen 1992). This 
minimum estimate of 23 bottleneck survivors is quite 
consetvative for a number of reasons. First, we used 
a conservative minimum estimate of the number extant in 
the early 1920s, based on counts of seals on land, whereas 
most seals tend to be at sea at any given time. Second, it is 
unlikely that the population grew at a sustained rate as high 
as 7% per year as the seals continued to be impacted by 
killing, harassment, and human occupation of terrestrial 
habitats into the 20th century (Ragen 1999). Thus, we think 
that considerably more than 23 seals survived in the early 
1890s. On average, effective population size is a tenth of the 
total population size (Frankham 1995), and Ne = 2 may 
likewise be underestimated. Finally, it is highly unlikely 
that the species remained at N. = 2 for 4 generations as 
there were at least 60 seals in 1893 and “numerous” seals by 
1915 (Le., less than 2 generations). Therefore, heterozygosity 


Table 3. Historic records of Hawaiian monk seal counts during or prior to the early 1920s (extracted from Hiruki and Ragen 1992) 


Location Year Count Comment 

French Frigate Shoals No data 

Laysan Island 923 4 2 of 4 were killed 
Lisianski Island 923 10 No census, 10 seal killed 
Pearl and Hermes Reef 923 125 Not full count of atoll 
Midway Atoll 913 5 ot 6 No exact count 

Kure 920 40 to 50 No exact count 

Total minimum 184 


The minimum count from each site was tallied to calculate minimum abundance during this postbottleneck period. 
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was likely less than 9% even prior to the 19th century 
bottleneck. 

Extremely low genetic diversity may be due to earlier 
hunting by indigenous Hawaiians and characteristics in- 
trinsic to the species. Average polymorphism and hetero- 
zygosity are lower in vertebrates than invertebrates, with 
mammals (especially carnivores) exhibiting the lowest 
genetic diversity (Nevo 1978). There is a pattern of reduced 
genetic variability with increased body size in mammals 
(Wooten and Smith 1985). Island populations may have 
lower genetic diversity than their mainland counterparts as 
a result of few founding individuals and because small 
population size increases the rate of genetic drift and 
decreases the accumulation of new mutations. For example, 
the Barrow Island population of the black-footed rock 
wallaby (Petrogale lateralis) has extremely low genetic diversity 
(H. = 0.053) as compared with mainland populations 
(H. = 0.56, 0.62), which is attributed to inbreeding over 
several generations (Eldridge et al. 1999). Contemporary 
inbreeding seems an unlikely cause of low genetic diversity 
in the Hawaiian monk seal because the heterozygosity deficit 
is low (Fig = 0.018). Likewise, we find no evidence for 
population structure (K = 1), which has been shown to 
deplete genetic diversity in small populations (Amos and 
Harwood 1998). 

Other large or carnivorous mammals have similarly 
undergone recent bottlenecks and have low levels of genetic 
diversity (Table 2), including those with total census sizes of 
less than 100 (northern hairy-nosed wombat, Lasiorhinus 
krefjtii) and 300 (North Atlantic right whale, Eubalaena 
glaciali) individuals. In such circumstances, some species 
have reduced fitness: wild and captive cheetah (Acinonyx 
jubatus) populations are often cited as examples of 
inbreeding depression. Genetic homogeneity in the cheetah, 
as demonstrated by the acceptance of skin grafts among 
unrelated individuals, may be correlated to high juvenile 
mortality, male fertility problems, and viral vulnerability 
(O’Brien et al. 1985, but see Merola 1994 and O’Brien 
1994). Relative to the cheetah (4 = 3.4; He = 0.39; 
Menotti-Raymond and O’Brien 1995), estimates of micro- 
satellite diversity in the Hawaiian monk seal (4 = 1.1; H. = 
0.026) raise the possibility of a genetic component to low 
juvenile survival and low reproductive rates. 

Perhaps the most appropriate comparisons of genetic 
diversity are within the family Phocidae, the true seals. 
Unfortunately, the Caribbean monk seal (Monachus tropicalis) 
went extinct in the mid-20th century (Kenyon 1977; 
LeBoeuf et al. 1986) before genetic diversity could be 
assessed. The Mediterranean monk seal (Monachus monachus) 
is the most threatened of pinniped species, with fewer than 
400 individuals remaining as a result of disease outbreak 
(Harwood 1998) and hunting dating back to ancient Greek, 
Roman, and Byzantine civilizations (Pastor et al. 2004). 
A single mitochondrial control region haplotype is shared by 
18 Mediterranean monk seals (Stanley and Harwood 1997). 
Microsatellite diversity is low but remains higher than that 
of the Hawaiian monk seal (Gemmell et al. 1997). Pastor 
et al. (2004) find 15 polymorphic loci of 39 tested (4 = 1.5; 
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H. = 0.16). Notably, the Hawaiian monk seal also has lower 
allelic diversity and heterozygosity than the northern 
elephant seal (M. angustirostris), the previous record for 
depleted genetic diversity among animals. 

Low genetic diversity has not prevented the northern 
elephant seal population from rebounding to more than 
175 000 (Weber et al. 2000), but monk seals have not been 
as resilient. Despite protection of the Hawaiian monk seal 
(through the US Endangered Species Act and the US 
Marine Mammal Protection Act) and an isolated, primary 
habitat in the Papahanaumokuakea Marine National Mon- 
ument, the species continues to decline at 4% per year 
(Carretta et al. 2007). The Mediterranean monk seal has 
endured centuries of exploitation, a recent disease-induced 
bottleneck (Harwood 1998) and smaller population sizes, yet 
maintains higher levels of genetic diversity, possibly a result 
of the preservation of alternate alleles by fragmentation of 
a once continuous range (Pastor et al. 2007). 

The Hawaiian monk seal sets a new standard for low 
genetic diversity in endangered species, enhanced, but 
apparently not solely caused, by a bottleneck in the 19th 
century. Since that time, monk seals have experienced 
robust population growth, thereby indicating that the 
species is viable despite minimal genetic variation. Recently, 
however, the species has declined as a result of low juvenile 
survival attributed to food limitation, shark predation, and 
entanglement in marine debris (Antonelis et al. 2006). 
Although genetic factors may not be driving the current 
trend, we cannot ignore theit potential impact on future 
population persistence because species with higher genetic 
diversity have experienced compromised fertility, reduced 
reproductive rates, high juvenile mortality, and disease 
epidemics. Despite these hazards, monk seals have stirred 
the waters of the Hawaiian archipelago for more than 10 
million years; given the opportunity, they will continue this 
ancient relationship. 
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